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Absence of zero-temperature transmission rate of a double-chain tight-binding model 
for DNA with random sequence of nucleotides in thermodynamic limit 
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The zero-temperature transmission rate spectrum of a double-chain tight-binding model for real 
DNA is calculated. It is shown that a band of extended-like states exists only for finite chain length 
with strong inter-chain coupling. While the whole spectrum tends to zero in thermodynamic limit, 
regardless of the strength of inter-chain coupling. It is also shown that a more faithful model for real 
DNA with periodic sugar-phosphate chains in backbone structures can be mapped into the above 
simple double-chain tight-binding model. Combined with above results, the transmission rate of 
real DNA with long random sequence of nucleotides is expected to be poor. 
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DNA is the biological molecule which keeps and propa- 
gates the secretes of life for all the creatures on the earth, 
and it has many interesting properties^. It consists of two 
chains of nucleotides. The nucleotides are of four types, 
usually denoted as A, T, C and G. Electrons are occu- 
pied on the nucleotides and can hop from one nucleotide 
to its intra-chain and inter-chain neighbors. Each pair of 
inter-chain nearest-neighbor nucleotides can only belong 
to one of the four 'base-pairs' A — T, T — A, C — G and 
G — C while other kinds of combinations are forbidden. 
Therefore, when the sequence of nucleotides in one chain 
of a DNA is determined, the sequence in the other chain 
is determined, too. Fig^a) shows schematically a part 
of a DNA sequence. 

The electronic transportation properties of DNA has 
attracted much attention^ since Elley et. al. suggested 
that DNA may become one-dimensional conductor 2 . 
Different conclusions are obtained by different exper- 
iments. Some experiments suggested that DNAs are 
conductors- Later, other experiment groups claimed 
that DNAs are insulators 4 . Direct measurement on single 
molecule showed that they seem to be semiconductors 5 . 
Some hints of superconduction behavior have even been 
reported^. Further experiments show that the conduc- 
tion behavior of DNA molecules of different base-pair 
sequences, say, identical base-pair sequence* and disor- 
dered base-pair sequence^, seems different. 

In the aspect of theoretical studies, a simple one- 
electron model was originally proposed by IguchiS. In 
this simple model, each nucleotide is represented by a site 
and different types of nucleotides have different on-site 
energies. Since the coupling between nucleotides is short- 
range, it is assumed that electron hopping occurs only be- 
tween nearest-neighbor nucleotides and the possible ran- 
domness in hopping constants is neglected. Thus the 
electron behavior in a DNA is represented by a double- 
chain tight-binding model(TBM). Many studies on elec- 
tronic properties of DNA are based on this model or its 
modified versions 1 ^^* 1 ^ 1 ^ 1 ^* 1 ^ 1 ^* 1 ^. In this short paper, 
we shall show that this simple model may not be appro- 



priate for real DNAs if they do have conductor behavior. 
A system should have extended states to provide conduc- 
tor behavior. It is well known by the famous Bloch the- 
orem that a periodic TBM has bands of extended states. 
However, the sequence of base-pairs in each chain of real 
DNAs are disordered rather than periodic, and this leads 
to random on-site energies in the corresponding TBM. 
According to the scaling theory of localization^!, uncorre- 
cted random on-site disorders in a one-dimensional TBM 
will drive all electronic states to exponentially localized 
states. At a first glance, the double-chain TBM for DNAs 
seems different from a common TBM with uncorrelated 
random on-site energies since the restriction of the com- 
bination of inter-chain base-pairs in DNA will induce lo- 
cal correlations between the two chains. However, by 
calculating the zero-temperature transmission rate spec- 
trum directly we shall show that this double-chain TBM 
for real DNAs cannot have conductor behavior, either. 
There have already been some theoretical studies on sim- 
ilar problemsiii 2 * 1 ^* 1 ^ 1 ^* 1 ^, and we shall also discuss and 
compare our results with them. 
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FIG. 1: (a)Schematic illustration of a part of DNA sequence; 
(b)the corresponding double-chain TBM of the sequence in 
(a). €a, er, ec, eo denote the different on-site energies for 
the corresponding nucleotides, t denotes intra-chain coupling 
and tinter inter-chain coupling. 



The Hamiltonian of the double-chain TBM for disor- 
dered DNAs is as following 
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e n! i and e ni 2 are the on-site energies of nucleotides in 
the two chains which can take four possible values e^, 
eT,£Cj e G- t and Unter are the intra-chain and inter-chain 
couplings, respectively. Fig^b) shows the correspond- 
ing double-chain TBM for the part of DNA sequence in 
FigUta). 

In order to introduce the effect of disorder in on-site en- 
ergies, we consider the simplest case that only two types 
of nucleotides, say, A and T, exist in the two chains while 
the sequence of A and T in each chain is random. This 
may be considered as a minimum model to include on-site 
disorder effect because in a general sequence in DNA all 
four possible base-pairs exist randomly which will only 
enhance the effect of disorders. For this special case, 
without loss of generality we can take the average on-site 
energy of nucleotides A and T as the energy reference 
point i.e., (ca + £t)/2 = 0, and set the intra-chain cou- 
pling t as the energy unit, i.e., t = 1. Then, we need 
only consider the case that 6a — +e and ex = — e, and 
e can be considered as the disorder strength in on-site 
energies. Fig|2Ja) shows a part of such a DNA sequence, 
and Fig[2Jb) is the corresponding TBM. 
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FIG. 2: The double-chain random tight-binding model for 
DNA. (a)Schematic illustration of a part of the double-strand 
DNA with only A and T as we consider; (b)the corresponding 
double-chain tight-binding model. 

The zero-temperature transmission rate of an eigen- 
state of energy E in unit of Gq = 2e 2 /h is given byM 
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cosh X n (E) 
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where X n (E) are eig envalues of the operator T^(E)f{E) 
and T(E) is the transfer matrix of the system. The eigen- 
values \ n (E) and the transmission rate spectrum can be 
calculated by standard transfer matrix algorithmic. This 
mathematical framework is essentially the same as the 
one Roche et. al^ have used in their study. We shall 
show that the transmission rate spectrum of the double- 
chain TBM for disordered DNAs is strongly suppressed 



and all electronic states are localized in the thermody- 
namic limit. 
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FIG. 3: The transmission rate spectrum g(E) for a double- 
chain tight-binding model with random on-site energies. 
(a)ti„ter = 1 with L = 10000, (b)rji, Iter = 3 with L — 
10000, (c)t lnter = 5 with L = 10000, {d)U nter = 9 with 
L = 10000, (e)t inter = 9 with L = 50000, {f)t inter = 9 with 
L = 100000. 

Let us look at the numerical results. Since all results 
are symmetric with E=0, we only provide the part of 
E > 0. Figgis for L = 10000 with (a)i mter = 1, (b) 
Unter = 3, (c) t mter = 5, (d) t inter = 9. The parameter 
for the disorder strength of on-site energies is taken as 
e = l, and the probabilities of the occurrence of both 
A and T in each chain are set as 1/2. Essentially the 
same results are obtained for other non-zero value of e 
and non-zero occurrence probabilities of A and T. 

For weak inter-chain coupling U n t er = 1, the transmis- 
sion rate is zero everywhere, which means that all states 
are localized. When inter-chain coupling increases, say, 
Unter = 3, peaks begin to exist in the spectrum, simi- 
lar to the behavior of a random single-chain TBM of fi- 
nite length. For further strong inter-chain coupling, say, 
Unter = 5, the curve begins to form a band with sharp 
peaks. With further increase of ti nter , say, ti nter = 9, 
a band of finite transmission rate emerges. However, 
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as shown in Fig[3Jd), (e) and (f) for ti nter = 9 with 
L = 10000, 50000 and 100000, respectively, it is clear- 
that this band tends to vanish at large L. Thus the exis- 
tence of a band of finite transmission rate in the case of 
strong inter-chain coupling does not mean the existence 
of truly extended states. It is only a finite-size effect 
which is expected to disappear when the chain length L 
tends to infinity. 

According to the above results, we may come to the 
following conclusions. A disordered double-chain TBM 
for real DNAs with weak inter-chain couplings does not 
have conducting electronic states. When inter-chain cou- 
plings become strong, the model with finite chain length 
forms an energy band of finite transmission rate, similar 
to that of a periodic single-chain TBM. However, this is 
due to finite-size effect and does not mean the existence of 
true conducting states. In the thermodynamic limit, the 
band tends to disappear and all electronic states are lo- 
calized. Therefore, the conductor behavior of real DNAs, 
if exists, cannot be explained by the simple one-electron 
double-chain TBM with random on-site energies. 
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FIG. 4: (a)Schematic illustration of the faithful model of a 
part of double-chain DNA with A and T nucleotides only; 
(b)the corresponding double-chain TBM of the sequence in 
(a). S and P denote sugar and phosphate sites, es and ep 
denote on-site energies of sugar and phosphate sites, respec- 
tively. €a and £t denote on-site energies of A and T nu- 
cleotides, respectively, tsp denotes intra-chain coupling be- 
tween a sugar and its neighboring phosphate. tpA and tpx the 
coupling between a phosphate and an A and T, respectively. 

However, it is well-known that there is a periodic 
sugar-phosphate chain in each backbone structure of real 
DNA. Therefore, a faithful model for a real double-chain 
DNA should be as shown in FigQJa) and (b), where S 
and P denote sugar and phosphate and A, T denote the 
nucleotides. Thus there comes a question on whether the 
simple model in Fig^ without sugar-phosphate periodic 
chains can be used for real DNA. Fortunately, we shall 
show below that the model in Fig0] can be renormal- 
ized into the model in Fig[3by the real-space decimation 
renormalization technique^. We shall still restrict to the 
simplest case that only A and T exist in the two chains. 

Let us first look at a single-chain DNA with series 
of random nucleotides as shown in Fig[S{a). We shall 
show that it can be mapped into a single chain Anderson 
TBM with random on-site energies. As a first step, let 
us renormalize a local phosphate-sugar-phosphate struc- 
ture into a phosphate-phosphate structure by decimation 
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FIG. 5: (a)Schematic illustration of the faithful model of a 
part of a single-chain DNA with only A and T nucleotides; 
(b)the corresponding single-chain TBM of the sequence in (a). 
S and P denote sugar and phosphate sites, es and ep denote 
on-site energies of sugar and phosphate sites, respectively, ea 
and er denote on-site energies for A and T nucleotides, re- 
spectively, tsp denotes intra-chain coupling between a sugar 
and its neighboring phosphate. tpA and tpp the coupling 
between a phosphate and A and T nucleotides, respectively. 



of the sugar site, as shown in FiglHl Perform this dec- 
imation process on all local phosphate-sugar-phosphate 
structures, one can decimate all sugar sites. The decima- 
tion process renormalizes the on-site energy of all phos- 
phate sites into e' P 



ep 



2t SP 



E 



(5) 



and introduces a direct coupling tpp between nearest- 
neighboring phosphates 



pp 



SP 



E-es 



(6) 



E is the eigen-energy we consider, es and ep are on- 
site energies of sugar and phosphate sites, and tsp is the 
coupling between neighboring sugar and phosphate sites 
in the original model of Fig[5] (Both e' P and tpp diverge 
at E = es which should be treated by considering directly 
the original stationary Schrodinger equations.) For each 
given E, the above decimation only gives a global shift 
for the on-site energies of each phosphate site and does 
not introduce any disorder because es, ep and tsp are 
constant values in the original model. 
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FIG. 6: A local phosphate-sugar-phosphate structure (a) is 
renormalized into a phosphate-phosphate structure (b) by 
decimation of the sugar site S. 

The second step is to renormalize a local phosphate- 
phosphate-phosphate structure by decimation of the nu- 
cleotide coupled with the middle phosphate as shown in 
Fig[3 This decimation process only renormalizes the on- 
site energy of the middle phosphate ep into 



e P = ep 



PA 



(7) 
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where tA is the on-site energy of the nucleotide coupled 
with the middle phosphate. (e' P diverges at E = ca, 
which means that the electron wave of energy E = ca 
is totally reflected at this point, similar to the 'anti- 
resonance' phenomenon^!) . Since the series of A and T 
nucleotide are random for the case we consider, this dec- 
imation process introduces effectively randomness into 
renormalized on-site energies of phosphate sites. Perform 
this decimation step for all local phosphate-phosphatc- 
phosphate structures, one finally obtain a single chain of 
phosphate sites with uncorrelated random on-site ener- 
gies. According to the scaling theoryii, the on-site dis- 
order localizes states of all eigen-energies. 
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FIG. 7: A local phosphate-phosphate-phosphate structure 
with the middle phosphate coupled with a nucleotide (a) is 
renomalized into a phosphate-phosphate-phosphate structure 
structure (b) by decimation of the nucleotide coupled with 
the middle P. 
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FIG. 8: A local structure (a) is renormalized into a local 
structure (b) by decimation of the nucleotide pair. 

Now let us consider the double-chain model as shown 
in FigQJb). Perform the same process as shown in FigO 
one can decimate all sugar sites. In order to map into 
the model of Fig[3 one needs to decimate all the nu- 
cleotide pairs. This can be done by considering the dec- 
imation process shown in Fig [5] where a pair of coupled 
nucleotides connecting two phosphates in each chain are 
decimated. This will renormalizc the on-site energies of 
the two phosphates as 
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(e'p is as in Eq.JSJ) and introduce a direct inter-chain 
coupling between the two phosphates as 



/ _ tpAtpTtAT 

tAT ~ (E-e A )(E-e T )-t A 



(9) 



(e' A , e' T and t' AT diverge at (E — €a){E — ex) — t AT = 
which should be treated by considering directly the orig- 
inal stationary Schrodinger equations.) It is easy to see 
that this decimation process introduces correlated ran- 
domness into on-site energies of each pair of inter-chain- 
coupled phosphate sites. Therefore, the model in FigQ] 
is mapped into essentially the same model as shown in 
Fig|21 According to numerical results for the model in 
Fig|21 we may conclude that no extended states exist in 
thermodynamic limit for both models. 
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FIG. 9: A double-chain TBM model for DNA with couplings 
between nucleotide and its neighboring sugar sites (a) is renor- 
malized into a double-chain TBM (b) by decimation of the 
sugar sites and nucleotide pairs. 

It should be noted that extended states may emerge 
when other kinds of couplings are included in the above 
model. For example, let us take into account the cou- 
plings between a nucleotide and its nearest-neighboring 
sugar sites, as shown in Fig[|Ja). Then, by decimation 
of both the sugar sites and nucleotide pairs, one can ob- 
tain a model of Fig^b) where both random on-site ener- 
gies and random couplings in each group of four nearest- 
neighboring sites are correlated. Such a model may have 
extended states because it looks similar to generalized 
dimer models which have been shown to have extended 
states in recent studies^. 

Before making the summary, we would like to 
make some discussions about our results and recent 
theoreticatiiii2ii2iiii£ii& and experimentafi&results on 
similar problems. Iguchii& has studied a single chain of 
DNA with periodic nucleotide series, which always has 
extended states due to Bloch's theorem. Roch et. ali^ 
have studied both periodic approximations of aperiodic 
DNA series and series extracted from real DNAs. Their 
results for series of real DNAs are similar with our re- 
sults. Yamadaii*i2ii£ and Yamada et. al>i^ have studied 
the influence of correlated on-site disorder in the model 
and found similar localization behaviors as our results 
suggest. In the aspect of experiments, Yoo et. ali have 
found that the conduction behavior of DNA with iden- 
tical base-pairs may be well explained by bands of ex- 
tended states separated by localized states which is pos- 
sible for a periodic TBM model. While Tran et. al.'s 
studies- on A— DNA with disordered sequence of base- 
pairs seem to suggest that no truly extended states exist, 
which agrees with our results. 

In summary, we have studied the disordered double- 
chain TBM with and without periodic sugar-phosphate 
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chain in the backbone structure for real DNAs. Numer- 
ical results combined with real-space decimation renor- 
malization technique show that in both cases no truly 
extended states exist in the thermodynamic limit, regard- 
less of how strong the inter-chain couplings are. There- 
fore, the double-chain TBM for real DNAs considered in 



this paper always has zero transmission rate in thermo- 
dynamic limit, and bad conduction is expected for real 
DNAs with long random series of nucleotides. 
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